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Notation 

a  -  Crack  Depth 

b  -  Beam  Width 

CMOD  -  Crack  Mouth  Opening  Displacement 

CTOD  -  Crack  Tip  Opening  Displacement 

E  -  Modulus  of  Elasticity 

G  -  Energy  Release  Rate 

KI  -  Stress  Intensity  Factor  for  Opening  Mode 

L  -  Beam  Length 

M  -  Moment 

P  -  Load 

P-Max  -  Maximum  Load 

R  -  Rotation  Factor 

w  -  Beam  Depth 

v  -  Poisson's  Ratio 


Chapter  1 
INTRODUCTION 

The  applicability  of  fracture  mechanics  to  concrete  has  been  under 
investigation  for  many  years.  Various  test  specimens  and  testing 
procedures  have  been  used  but  few  have  had  successful  and  consistant 
results.  Concrete  has  proved  to  be  a  notch  sensitive  material,  that  is, 
it  behaves  differently  when  notched  with  teflon  or  a  sawcut,  than  it  does 
when  it  is  precracked.  This  factor  alone  affects  many  of  the  previous 
investigations  since  many  were  done  with  notched  beams.  Many  methods  of 
analysis  have  been  proposed,  as  have  many  testing  procedures.  In  spite  of 
the  problems  encountered  it  is  still  felt  that  fracture  mechanics  may  be 
the  best  way  to  describe  the  fracture  behavior  of  concrete. 

An  extensive  study  has  been  started  in  an  attempt  to  verify  the 
relationships  between  the  fracture  parameters  for  plain  concrete  in 
bending.  The  results  presented  here  are  representative  of  small  beams 
with  a  width  of  3  in.,  a  depth  of  4  in.,  and  a  span  of  15  in.  loaded  in 
three-point  bending.  A  modified  compliance  calibration  technique  was  used 
to  precrack  the  beams  and  load  versus  crack  mouth  opening  displacement 
(CMOD)  and  load  versus  load  point  displacement  (LPD)  plots  were  obtained 
simultaneously.  Five  precracked  beams  of  each  crack  depth  to  beam  depth, 
a/w,  ratio  of  0.3,  0.5,  and  0.7  were  failed  as  well  as  two  beams  each  of 
the  same  a/w  ratios  with  teflon.  A  detailed  description  of  the  methods 
used  is  found  in  Chapter  3. 


The  results  obtained  included  stress  intensity  factors  (Kjr),  energy 
release  rates  based  on  KjCj  energy  release  rates  based  on  the  J-integral 
concept,  energy  release  rates  found  from  an  energy  method  adapted  from 
Petersson  (4),  and  the  crack  tip  opening  displacement  (CTOO).  The  use  of 
the  crack  length  at  the  point  of  instability  was  compared  to  the  use  of 
the  precracked  length  as  determined  by  the  dye.  These  results  are 
presented  in  Chapter  4. 

Results  obtained  show  fairly  consistant  energy  release  rates  based  on 
the  energy  approach  and  the  J-integral  approach,  especially  for  larger 
crack  lengths.  Those  based  on  the  stress  intensity  factors  tended  to  be 
much  lower  than  those  found  by  the  other  two  methods.  The  CTOD  values 
however,  show  consistant  results  for  the  extended  crack  lengths.  The 
conclusions  and  summary  are  found  in  Chapter  5. 


Chapter  2 
LITERATURE  REVIEW 

Many  experiments  have  been  performed  in  an  attempt  to  prove  that  the 
fracture  toughness  of  concrete  is  a  material  property.  Factors  such  as 
test  specimen,  testing  procedure  and  method  of  analysis  have  been 
investigated  and  some  of  the  pertinent  investigations  are  described 
below. 

An  intensive  study  by  Go  and  Swartz  (3)  investigated  four  criteria 
which  were  felt  to  be  related  to  concrete  failure.  They  were  energy 
release  rate,  J-integral  concept,  stress  intensity,  and  the  crack  tip 
opening  displacement  (CTOU).  Their  studies  also  included  the  feasibility 
of  the  compliance  calibration  technique  and  considered  crack  growth  and 
the  influence  of  the  microcracking  zone. 

Studies  on  the  stress  intensity  factor  consisted  of  a  comparison 
between  the  bending  analogy,  finite  element  method,  and  the  Srawley 
formula  in  determining  the  stress  intensity  factor.  An  equation  for 
estimating  Kj  was  derived  using  the  least  squares  method  and  is  as 
follows  (3): 

Kj  »— £-( Az^Bz+C+Dz^+Ez"2  ) (2.1) 

bw1-5 

M  =  moment  at  midspan 
b  •  beam  width 
w  =  beam  depth 
z  =  1  -  a/w 

A,  B,  C,  0,  E  =  Constants  determined  for  different 
span/depth  ratios 


The  applicability  of  the  compliance  calibration  method  as  it  relates 
to  concrete  was  investigated  by  comparing  crack  length  estimates  found 
using  the  compliance  curve  with  actual  crack  lengths  revealed  by  the  dye 
penetrant  (3).  Thirty  specimens,  made  in  two  sets,  were  three  inches 
wide,  four  inches  deep  and  had  a  fifteen  inch  span.  Two  beams  per  set 
were  notched  at  midspan  to  various  depths  and  used  to  determine  compliance 
and  a  plot  of  compliance  versus  a/w  obtained.  Fourteen  of  the  beams  had 
teflon  inserts  and  were  loaded  to  failure  under  load  control.  A  load 
versus  CMOD  trace  was  obtained  and  the  compliance  found  from  the  initial 
slope.  Estimated  crack  lengths  from  the  compliance  curve  were  compared 
with  the  actual  crack  lengths.  The  remaining  twelve  specimens  were 
notched  0.5  in.  (12.7  mm)  then  loaded  under  strain  control  to  generate  a 
crack  of  desired  length  using  the  compliance  curve.  Dye  was  inserted  and 
worked  into  the  crack  while  the  beam  was  in  an  inverted  position.  The 
beam  was  then  righted  and  loaded  to  failure  under  load  control.  The 
estimated  crack  length  was  determined  using  the  compliance  curve  and  the 
slope  of  the  trace  after  crack  closure  was  overcome.  This  crack  depth  was 
compared  with  the  actual  crack  depth  revealed  by  dye. 

Go  and  Swartz  (3)  concluded  that  application  of  the  compliance 
calibration  technique  should  be  associated  with  the  real  crack  length  of  a 
precracked  beam  instead  of  a  sawcut  depth.  To  determine  energy  release 
rates,  sixteen  specimens  were  tested  in  three-point-bending  and  a  trace  of 
either  load  versus  CMOD  or  load  versus  LPD  and  a  compliance  curve  based  on 
LPD  was  obtained. 

The  stress  intensity  approach  considered  the  specimen  to  be  in  plane 
strain  and  assumed  that  the  sliding  mode  and  tearing  mode  stress  intensity 


factors  were  negligible  compared  to  the  opening  mode.  The  equation  for 

plane  strain  then  reduced  to: 

l-v2 
s--£-Kt2 (2.2) 

in  which  Kj  was  based  on  the  extended  crack  length. 


The  J-integral  approach  is  based  on  the  concept  that  for  an  infinitesimal 
amount  of  crack  extension,  the  decrease  in  stored  elastic  energy  of  a 
cracked  body  loaded  under  displacement  control  is  identical  to  the 
decrease  in  potential  energy  when  loaded  under  load  control.  To  use  this 
concept  experimentally  ,  a  plot  of  the  energy  required  to  trigger  the 
instability  versus  the  unextended  a/w  was  made  and  the  slope  of  this 
straight  line  gave  the  energy  change  per  change  in  a/w.  This  energy  rate 
was  then  divided  by  the  effective  remaining  area,  that  is,  a  surface  15£ 
larger  than  the  area  revealed  by  dye  to  account  for  the  roughness  of  the 
surface  (3).  A  modification  of  Petersson's  (4)  approach  was  also 
investigated.  The  total  energy  consumed  was  calculated  by  finding  the 
area  up  to  the  point  of  instability  of  the  failure  plot.  The  energy 
release  rate  was  then  calculated  by  dividing  the  energy  by  the  effective 
remaining  area  at  the  point  of  unstable  crack  growth. 

Conclusions  reached  by  Go  and  Swartz  were  that  the  average  energy 
release  rate  for  the  J-integral  and  stress  intensity  approaches  based  on 
stable  crack  growth  and  microcracking  extension  were  in  good  agreement. 
Petersson's  method,  when  not  considering  stable  crack  growth,  gave  lower 
values.  However,  good  agreement  with  the  other  two  methods  was  found  when 
Petersson's  method  was  applied  considering  crack  extension.   The  results 


based  on  the  assumption  of  zero  crack  growth  are  clearly  a  lower  bound  on 
the  actual  failure  results  either  in  terms  of  energy  release  rate  or 
stress  intensity  according  to  Go  and  Swartz.  It  was  also  determined  that 
linear-elastic  fracture  mechanics  concepts  are  valid  in  the  study  of 
concrete  fracture. 

The  COD  approach  for  specifying  fracture  toughness  was  examined. 
Based  on  the  concept  that  the  rotating  center  of  the  crack  is  at  the 
strain  reversal  point  (i.e.  neutral  axis),  Go  and  Swartz  used  William's 
stress  function  to  evaluate  the  rotation  factor  R  for  the  3  point  bending 
specimen  (3).  The  rotation  factor,  R,  is  defined  as  shown  in  Figure  2.1. 
For  span/depth  ratios  equal  to  3.75,  the  average  value  of  R  was  0.45  with 
maximum  deviation  less  than  5%.  Using  this  rotation  factor,  the 
relationship  between  the  CMOD  and  CTOD  was  found  to  be 

S  (CTOU)  ■   R(w-a)   V(CMOD) (2.3) 

R(w-a)+a+z 

for  beams  in  three  point  bending.  The  constant,  z,  is  the  thickness  of 
the  knife  edge  as  shown  in  Figure  2.1.  It  was  concluded  that  the  CTOU 
showed  promise  as  a  fracture  criterion  for  concrete.  It  may  be  a  valid 
alternative  treatment  instead  of  G,  or  Kt. 

A  study  done  by  Fartash  (2)  was  conducted  to  compare  the  behavior  of 
beams  with  notches  to  those  with  natural  cracks.  A  total  of  ninety  six 
beams  with  a  width  of  3  in.,  a  depth  of  4  in.  and  a  span  of  15  in.  were 
made  and  tested  in  3  and  4  point  bending.  Twelve  beams  were  used  for 
compliance  curves,  forty-eight  had  teflon  inserts,  and  forty-eight  were 
precracked.  The  crack  lengths  used  in  this  study  gave  a/w  ratios  of  0.3, 


0.5,  0.7.  Of  the  specimens  used  in  this  investigation  half  of  them  had  a 
strength  of  3200  psi  and  half  had  a  compressive  strength  of  6700  psi. 

All  testing  was  done  using  an  electro-hydraulic  materials  testing 
system  (MTS).  Using  the  compliance  technique  developed  by  Swartz,  Hu  and 
Jones  (5),  the  plain  beams  were  initially  notched  .40  in.  (10.16  mm)  to 
.45  in.  (11.43  mm)  then  precracked  to  the  desired  length.  This  was  done 
by  loading  the  beam  in  strain  control  until  the  slope  of  the  load 
displacement  plot  decreased  at  which  point  the  load  was  removed.  The 
compliance  of  the  straight  part  of  the  curve  was  used  to  determine  the 
estimated  crack  length  from  the  compliance  calibration  curve.  This  was 
repeated  until  the  desired  crack  length  was  reached.  All  precracked 
specimens  were  then  loaded  to  failure  using  load  control.  The  teflon 
beams  were  loaded  to  0.6  P-max  to  determine  the  compliance  based  on  the 
CMOD  trace,  then  unloaded  and  reloaded  to  failure. 

A  load  versus  CMOD  was  obtained  for  each  specimen  tested.  From 
these,  compliance  values  were  calculated  as  well  as  net  bending  stresses. 
The  average  bending  stress  for  each  a/w  ratio  of  approx.  0.3,  0.5,  0.7  was 
plotted  versus  the  a/w  ratio.  The  plots  of  average  stress  intensity 
values  versus  average  a/w  ratios  indicated  that  both  batches  yielded 
notch-sensitive  specimens.  It  was  found  that  the  value  of  stress  intensity 
decreased  with  increasing  a/w  ratios.  A  relationship  between  stress 
intensity  values  for  precracked  and  prenotched  beams  for  various  a/w  ratios 
was  found  to  be: 


—  ■  .9137+1. 6118 (a/w)  (2.4) 

KIT 


Fartash  concluded  that  in  all  cases  the  naturally  cracked  beams  yield 
higher  failure  loads,  stress  intensity  values,  and  bending  stresses  than 
notched  beams  with  the  same  crack  length.  True  fracture  toughness  cannot 
be  obtained  with  notched  beams  unless  a  correlation  between  precracked  and 
notched  beams  is  used. 

The  presence  of  the  fracture  zone  in  concrete  was  investigated  by 
Petersson  (4).  Stable  tensile  tests  were  performed  on  concrete  specimens 
with  a  maximum  aggregate  size  of  5/16  in.  (8  mm.)  and  a  water  cement 
ratio  of  0.6.  After  being  cast  in  steel  molds,  the  specimens  were  cured 
in  lime  saturated  water  until  one  day  before  testing  when  they  were 
removed  and  wrapped  in  plastic  foil.  Total  curing  time  was  seven  days. 
The  specimens  were  allowed  to  dry  slightly  while  being  placed  in  the 
testing  machine,  then  all  sides  but  one  were  rewetted.  As  the  fracture 
zone  developed  it  became  visible  as  a  dark  area  on  the  dry  side  of  the 
speciman.  This  dark  area  was  a  result  of  the  water  absorbed  by  the  zone 
due  to  capillary  forces. 

Petersson  observed  that  the  fracture  zone  was  about  9.84  X  10-4  in. 
(25  X  10-3  m)  to  19.69  X  10-4  in.  (50  X  10-3  ran)  wide  when  the  crack 
became  visible.  The  width  of  the  fracture  zone  appeared  to  be  of  the  same 
magnitude  as  the  size  of  the  maximum  aggregate  particle. 

In  addition  to  investigating  the  fracture  zone  Petersson  also 
determined  the  fracture  energy  (Gp)  which  he  defined  t0  be  the  mount  of 
energy  necessary  to  create  one  unit  of  area  of  a  crack.   In  order  to 


measure  this  accurately,  the  speciman  was  loaded  to  failure  under 
displacement  control  to  obtain  a  stable  and  well  defined  length  of  crack 
propagation.  Notched  beams  were  loaded  in  3-pt-bending  and  a  stable  load 
to  failure  plot  obtained.  In  half  of  the  beams  weights  were  placed  on  the 
beam  outside  the  supports  to  compensate  for  the  energy  due  to  the  weight 
of  the  beam  alone.  These  results  were  compared  to  those  obtained  without 
compensation  for  the  weight  of  the  beam  and  it  was  found  that  the  true 
value  obtained  by  compensation  gave  values  50%  to  250%  higher  than  the 
uncompensated  beams.  Two  methods  of  estimating  the  true  Gp  hased  on  the 
results  of  uncompensated  beams  were  evaluated.  Additional  testing  was 
done  on  prenotched  (teflon  inset)  and  notched  (sawcut)  beams  as  well  as 
direct  tension  specimens.  The  specimens  were  loaded  in  displacement 
control  and  the  plot  of  load  versus  displacement  obtained.  The  Gp  values 
were  calculated  as  the  area  under  the  stable  curves  divided  by  the 
uncracked  area  at  the  beginning  of  testing.  In  determining  the  energy 
consumed  by  the  fracture  zone,  the  amount  of  energy  consumed  by  the 
material  outside  the  fracture  zone  before  the  tensile  strength  was  reached 
was  subtracted.  This  was  done  by  drawing  a  straight  line,  parallel  to  the 
major  slope  of  the  increasing  portion  of  the  graph,  to  the  point  of 
maximum  load.  The  area  between  this  line  and  the  increasing  portion  of 
the  curve  was  the  area  subtracted. 

Petersson  concluded  that  although  Gp  is  slightly  affected  by  beam 
depth,  it  seems  to  be  useful  as  a  material  parameter.  However  the  Gp  must 
be  calculated  using  the  stable  plot  obtained  using  displacement  control 
since  dynamic  effects  are  incurred  when  using  load  control  causing  crack 
growth  to  become  unstable. 
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FIG.  2.1    RELATIONSHIP    EETwEE\|    CTOD    AND   CMOD 
FOR  A   RIGID  ROTATION    ABOUT  THE  CENTER    OF    ROTATION 
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Chapter  3 
EXPERIMENTAL  PROGRAM 

3.1  Test  Specimens 

One  size  of  beam  was  constructed  to  the  following  dimensions: 

1  =  16  in.  (406.4  mm) 
w  •  4  in.  (101.6  mm) 
b  =  3  in.  (76.2  ran) 

Figure  3.1  shows  the  beam  dimensions.  One  mix  design  was  used  as  given  in 
Table  3.1.  This  gave  a  nominal  cylinder  compressive  strength  of  3100  psi 
(55.81  MPa),  Poisson's  Ratio  of  .195,  and  a  Modulus  of  Elasticity  of  5.34 
X  105  psi  (36.79  X  103  MPa).  Figures  3.2  and  3.3  show  stress  versus 
transverse  strain  and  stress  versus  longitudinal  strain  respectively. 
Cylinder  data  are  listed  in  Appendix  IV  with  the  raw  data.  Two  sets  of 
twenty  beams  each  were  cast  for  a  total  of  forty  beams.  Care  was  taken  in 
construction  and  curing  of  the  beams  to  ensure  that  the  two  sets  would 
have  nearly  identical  material  properties  at  the  time  of  testing. 

Of  the  forty  beams  cast,  fourteen  were  used  for  compliance 
calibration:  two  beams  for  each  a/w  ratio  of  .3,  .4,  .5,  .6,  .7,  .8,  and 
.9.  Fifteen  beams,  five  of  each  a/w  ratio  of  .3,  .5,  and  .7,  were 
precracked  and  six  beams,  two  of  each  a/w  ratio  of  .3,  .5,  and  .7,  had  a 
teflon  insert  which  simulated  the  a  real  crack  without  the  effect  of 
aggregate  interlock.  The  .003  in.  (.076  mm)  width  teflon  strips  were 
inserted  using  the  method  described  in  Reference  2.  One  beam  was 
instrumented  with  strain  gages  in  an  attempt  to  determine  the  strain 
distribution  ahead  of  the  crack  tip.  The  remaining  beams  were  used  as 
"spares". 
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3.2  Testing  Machine  and  Set  Up 

An  electro-hydraulic  materials  testing  system  (MTS)  was  used 
throughout  the  testing  program.  Crack  mouth  opening  displacements  (CMOO) 
and  load  point  displacements  (LPD)  were  monitered  simultaneously  through 
the  use  of  commercially  available  displacement  transducers  (MTS  632.05 
B-60)  which  have  a  maximum  sensitivity  of  +  .002  in.  (  +  .0508  mm)  per  10 
volt  full  scale  output.  During  loading  of  the  specimens,  simultaneous 
traces  of  load  versus  CMOD  and  load  versus  LPD  were  obtained. 

The  MTS  made  it  possible  to  load  the  specimens  in  three  different 
control  modes.  Under  load  control,  the  span  responds  to  the  amount  of 
load  as  the  primary  feedback.  Operating  in  this  control  mode  assures  a 
constant  load  rate  regardless  of  the  rate  of  crack  propagation  and  was 
used  when  loading  the  beams  to  failure.  Under  displacement  control,  the 
span  responds  to  the  displacement  of  the  CMOD  transducer  as  the  primary 
feedback.  In  this  control  mode,  it  is  possible  to  crack  the  beam  to  a 
desired  depth  without  much  danger  of  premature  failure  since  the  rate  of 
CMOD  is  controlable  instead  of  the  load  rate.  Precracking  of  the  beams 
was  done  under  displacement  control.  The  third  control  mode  is  stroke 
control.  It  uses  the  displacement  of  the  loading  head  as  its  primary 
feedback.  This  mode  was  used  in  the  dye  application  procees  due  to  its 
sensitivity  and  controlability. 

Appropriate  scale  settings  on  the  plotters  were  important  to  obtain 
the  best  size  traces.  The  plotters  used  were  MTS  431. 13A  -  02  (Type  200 
Control  Module)  and  a  summary  of  X  (load)  and  Y  (displacement)  axis  scale 
settings  follows: 
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X-axis    Metric  Setting 

Ranges  using  Calib.  setting: 

0.5%  per  cm  =  0.05  V/cm 
1.0%  per  cm  =  0.10  V/cm 
2.5%  per  cm  =  0.25  V/cm 
5.0%  per  cm  =  0.50  V/cm 
10.0%  per  cm  =  1.00  V/cm 

CM00  -  Range  2  +  1  X  10-2  in./10V  =  +  1  X  10-3  tn./V 

0.5%  :  1  cm  =  5.0  X  10"5  in. 

1.0%  :  1  cm  »  1.0  X  10-4  tn. 

2.5%  :  1  cm  =  2.5  X  10"*  in. 

5.0%  :  1  cm  =  5.0  X  10-4  1n. 

CMOD  -  Range  1  _+  2  X  10-2  in./10V  =  +  2  X  10"3  in./V 

0.5%  :  1  can  ■  1  X  10-4  ,„. 

1.0%  :   1  cm  =  2  X  10-4  j„t 

2.5%  :   1  cm  =  5  X  10"4  in. 

5.0%  :   1  cm  =  1  X  10-3  in. 

LPO  -  Range  2  +  9.72  X  10-4  in./v 


0.5%  :  1  cm  =  4.86  X  10"5  in. 

1.0%  :  1  cm  =  9.72  X  10"5  in. 

2.5%  :  1  cm  ■  2.43  X  10-4  in> 

5.0%  :  1  cm  »  4.86  X  10"4  in. 


LPD  -  Range  1     +  19.08  X  10"4  in./V 


0.5%  :  1  cm  -  9.54  X  1Q-5  in. 

1.0%  :  1  cm  =  19.08  X  10-5  in. 

2.5%  :  1  cm  =  4.77  X  10"4  in. 

5.0%  :  1  cm  -  9.54  X  10-4  in. 


Y-axis    Metric  Setting 

Ranges  using  Calib.  setting: 

0.5%  per  cm  =  0.05  V/cm 

1.0%  per  cm  =  0.10  V/cm 

2.5%  per  cm  =  0.25  V/cm 

5.0%  per  cm  =  0.50  V/cm 

10.0%  per  cm  =  1.00  V/cm 

CMOD  and  LPO  -  All  Ranges  1.0  V  =  1000  lbs. 
Using  load  cell  with  X  10 

0.5%  :  1  cm  =  50  lb. 
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1.0*  :  1  cm  ■  100  lb. 
2.5%  :  1  cm  =  250  lb. 
5.0%  :  1  cm  =  500  lb. 


The  settings  used  are  listed  on  each  trace. 

The  rate  at  which  the  beams  are  loaded  is  governed  by  both  the  span 
and  the  frequency  controls  on  the  MTS.  The  span  was  generally  set  on  110 
and  the  frequency  at  1.0  with  the  function  generator  set  at  0.1/1.1  on  the 
MTS  model  panel  410.21.  During  loading  to  failure,  this  setting  gave  an 
average  loading  rate  of  19  lb. /sec.  (84.55  N/sec.) 

A  listing  of  the  settings  used  throughout  the  testing  is  as  follows: 

Model  410.21   Panel 

Function:  upward  ramp  function 
Function  Generator:  0.1/1.1 
Frequency:  1.0 

406  Controller 

Cal  factor  =  4.21 

Excitation  ■  4.41 

Gain  =  8 

Rate  =  4.6 

P  =  0 

Fdbk  Select  =  XDCR1  =  stroke  control 
XDCR2  =  load  control 
EXT  ■  strain  control 

Problems  within  the  MTS  caused  delays  in  testing  and  an  erratic 
testing  schedule. 

3.3  Testing  Procedure 

(a)  Compliance  Calibration 

Each  compliance  specimen  was  initially  notched  at  midspan  using  a 
concrete  saw.  The  notch  was  0.13  in.  (3.302  mm)  wide  and  eighty  percent 
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of  the  desired  crack  length  in  depth.  The  purpose  of  the  starter  notch 
was  to  ensure  cracking  at  midspan.  After  the  beam  had  been  notched,  the 
yoke  shown  in  Figure  3.4  and  the  CMOD  transducer  were  mounted  on  the 
specimen  and  placed  in  the  loading  set  up  as  shown  in  Figure  3.1.  The 
beam  was  then  loaded  in  three  point  bending  using  displacement  control  and 
a  plot  of  load  versus  CMOD  was  obtained.  The  inverse  slope  of  the  major 
portion  of  this  curve  gave  the  compliance  value  for  the  beam's  current  a/w 
ratio.  Loading  was  continued  until  the  slope  began  to  decrease  indicating 
crack  growth.  Load  was  removed  then  reapplied,  and  the  compliance  was 
measured  from  the  new  trace.  This  procedure  was  repeated  until  it 
appeared  that  the  desired  crack  length  had  been  reached.  The  CMOD 
transducer  was  removed  and  the  specimen  was  inverted  and  placed  in  a 
reverse  three  point  bending  configuration  as  shown  in  Figure  3.5.  Using 
stroke  control,  the  beam  was  loaded  to  approximately  seventy-five  percent 
of  the  last  precracking  load  applied.  This  opened  the  crack  for  dye 
insertion.  Vanish,  a  product  of  The  Drackett  Products  company,  was 
inserted  using  a  pipette  and  the  load  then  cycled  to  work  the  dye  into  the 
crack.  Once  the  dye  had  been  worked  in,  the  CMOD  transducer  was  attached 
again,  and  the  specimen  replaced  in  the  loading  apparatus  as  during 
precracking.  The  loading  head  was  raised  to  where  it  touched  the  upper 
head,  but  no  load  was  being  applied.  At  this  point,  the  LPD  transducer 
was  attached  to  the  system  shown  in  Figure  3.6.  The  beam  was  then  loaded 
to  failure  using  load  control  and  simultaneous  traces  of  load  versus  CMOD 
and  load  versus  LPD  were  obtained. 

The  actual  average  crack  depth  was  determined  by  finding  the  area  of 
the  beam  penetrated  by  the  dye  and  dividing  it  by  the  width  of  the  beam. 
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Typical  failure  surfaces  are  shown  in  Figure  3.7.  The  crack  depth  found 
corresponds  to  the  compliance  of  the  failure  curve.  This  data  provides 
one  point  on  the  compliance  curve. 

Each  of  the  fourteen  compliance  specimens  were  tested  in  this  manner 
and  a  compliance  versus  a/w  ratio  curve  drawn. 

(b)  Precracked  Beams 

Each  precracked  beam  was  initially  notched  with  a  sawcut  as  the 
compliance  specimens  were.  A  crack  was  then  generated  to  the  desired 
depth  by  loading  the  beam  until  the  corresponding  compliance  value  as 
found  on  the  compliance  curve,  was  reached.  As  previously  described  the 
dye  was  applied  and  the  beam  loaded  to  failure.  Traces  of  load  versus 
CMOD  and  load  versus  LPD  were  obtained  as  well  as  the  areas  of  the  failure 
surfaces.  Typical  traces  are  shown  in  Figure  3.8  and  3.9.  In  all  traces 
the  vertical  axis  gives  the  applied  load.  The  horizontal  axis  gives  LPD 
or  CMOD  as  appropriate.  The  scale  factors  are  shown  on  each  trace. 

(c)  Beams  With  Teflon  Insert 

Since  no  precracking  and  dye  insertion  were  needed,  the  beams  were 
instrumented  with  both  transducers  and  loaded  to  failure  under  load 
control.  Traces  of  load  versus  CMOD  and  load  versus  LPD  were  obtained. 
The  cracked  area  caused  by  the  teflon  was  readily  visible. 

(d)  Beam  With  Strain  Gages 

The  beam  was  initially  sawcut  0.96  in.  (24.384  mm)  to  ensure  cracking 
at  midspan.  Fifteen  gages,  EA-O6-120LZ-120  with  a  gage  length  of  .375  in. 
(9.53mm),  were  affixed  to  each  side  of  the  beam  at  midspan  above  the  crack 
tip.  Strain  data  was  obtained  and  recorded  through  the  use  of  an  Optim 
data  aquisition  system.  The  beam  was  loaded  using  displacement  control 
until  cracking  occured.  Strain  readings  were  taken  and  the  compliance  at 
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that  point  measured.   This  procedure  was  repeated  for  increased  crack 
depths.  Figure  3.10  shows  the  gage  placement. 
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FIG.  3.1      THREE     POINT    BENDING     CONFIGURATION 
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FIG.  3.5      REVERSE   THREE     POINT  BENDING    CONFIGURATION 
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FIG-  3.6      LPD     TRANSDUCER    SETUP 
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A/W  =  0.326 
BEAM   NO.  C-2 


A/W=  0.525 
BEAM    NO.  C-4 


A/w  =  0.685 

EEAM    NC.  B-19 


FIG,  3.7       TYPICAL     FAILURE    SURFACES 
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Table  3.1  Mix  Design 

Batch  1  Batch  2 

Water/Cement                 0.50  0.50 

Cement  Type                   I  I 

S.G.  Sand*                   2.65  2.65 

S.G.  Aggregate*               2.56  2.56 

S.G.  Cement                  3.15  3.15 

%  Sand  by  weight              32.68%  32.68% 

%  Aggregate  by  weight          47.46%  47.46% 

%  Cement  by  weight             13.24%  13.24% 

%  Water  by  weight              6.62%  6.62% 

Density  of  concrete      149.7  pcf  (23.48  kN/m3)  149.7  pcf  (23.48  kN/m3) 

Curing  Time                 145  days  138  days 

Compressive  Strength       7950  psi  (54.78  MPa)  8130  psi  (56.02  MPa) 

Tensile  Strength**         601  psi  (4.14  MPa)  665  psi  (4.58  MPa) 

Superplasticizer               400  ml  300  ml 

Slump                  7.25  in.  (184.15  mm)  7.00  in.  (177.80  mm) 

Sand  Fineness 

Modulus                   2.91  2.91 

Maximum  Aggregate         0.75  in.  (19.05  mm)  0.75  in.  (19.05  mm) 
Size 

• 
Sand  and  aggregate  properties  and  testing  procedures  in  Appendix  I. 
** 

Tensile  strength  determined  by  split  cylinder  test  as  described  in 
ASTM. 
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Chapter  4 
EXPERIMENTAL  RESULTS 

4.1  Compliance  Beams 

The  fourteen  compliance  beams  were  tested  as  described  previously. 
Two  compliance  curves  were  drawn  by  plotting  the  compliance  value  versus 
the  a/w  ratio.  One  was  based  on  compliance  values  found  from  the  load 
versus  CMOD  traces  and  the  other  was  based  on  values  found  from  the  load 
versus  LPD  traces.  These  compliance  values  were  calculated  using  the 
slope  of  the  major  portion  of  the  trace  rather  than  the  initial  slope 
since  the  initial  slope  showed  evidence  of  crack  closure  and  slippage  of 
the  testing  apparatus.  The  slippage  referred  to  was  primarily  in  the  LPD 
set  up  and  is  noticeable  in  many  of  the  graphs.  A  load  versus  LPD  trace 
without  slippage  and  one  with  slippage  are  shown  in  Figure  4.1. 
Compliance  curves  are  shown  in  Figures  4.2  and  4.3. 

Some  difficulty  was  had  in  precracking  the  compliance  specimens  with 
a/w  ratios  of  .7,  .8,  and  .9. 

4.2  Precracked  Beams 

The  fifteen  precracked  beams  were  tested  as  previously  described. 
Stress  intensity  values  and  energy  release  rates  were  calculated  for  both 
the  unextended  and  the  extended  crack  lengths.  The  value  of  the  crack 
tip  opening  displacement  (CTOD)  was  also  obtained. 

(a)  Determination  of  Extended  a/w 

The  extended  a/w,  corresponding  to  the  point  of  unstable  crack 
growth,  was  calculated  by  drawing  a  line  from  the  origin  of  the  CMOD  trace 
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to  the  point  at  which  the  load  begins  to  drop  off,  that  is,  the  point  of 
unstable  crack  growth.  The  compliance  based  on  the  slope  of  this  line  was 
calculated  and  the  extended  a/w  found  using  the  CMOD  compliance  curve. 
The  CMOD  compliance  curve  was  used  since  it  gave  more  consistant  results 
than  the  LPD  compliance  curve,  however  calculations  were  made  using  the 
LPO  values  so  a  comparison  could  be  made. 

(b)  Stress  Intensity  Values 

The  stress  intensity  values  were  calculated  using  the  formula 
developed  in  Reference  3  by  Go  and  Swartz.  A  least  squares  method  was 
used  to  estimate  Kj  with  errors  less  than  3%.  The  formula  is  shown 
below: 

Kj  =  _  M   (Az2  +  Bz  +  C  +  Dz"1  +  Ez"2) (2.1) 

bw1.5 

M  =  Midspan  bending  moment  at  failure 

b  ■  beam  width 

w  ■  beam  depth 

z  =  1  -  a/w 

A  =  -0.065* 

B  =  -3.483* 

C  =  -0.120 

D  =  5.706* 

E  =  0.166* 

*  Coefficients  for  L/w  =  3.75  as  found  in  Reference  3. 
Stress  intensity  values  were  calculated  for  unextended  values  of  a/w  and 
for  extended  values  of  a/w  (KIc).  These  values  are  listed  ln  Table  4-1> 
4.2,  and  4.3. 

Average  values  of  Kj  for  unextended  crack  length  were  8.41  x  102 
lb. -in. -3/2  (g.249  X  105  N-m-3/2),  7.66  x  102  lb. -in. -3/2  (8.417  X  lfl5 
N-m-3/2),  5.62  x  102  lb. -in. "3/2  (6.177  X  105  N-m-3/2)  for  approximate  a/w 
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ratios  of  0.3,  0.5,  0.7  respectively.  Average  values  of  <,„  for  extended 
crack  lengths  based  on  CMOD  were  12.07  x  102  lb.-in.~3/2  (13.27  X  105 
N-nT^2),  9.44  X  102  lb.-in.~3/2  (10.38  X  105  N-nf3/2),  8.78  x  102 
lb. -in."3''2  (9.65  X  105  N-m"3/'2),  for  the  same  respective  a/w  ritios. 
Average  values  of  KIC  for  extended  crack  lengths  based  on  LPD  were  10.30  X 
102  lb.-in."3/2  (11.43  X  105  N-nT3/2),  7.92  X  102  lb.-in._3/2  (8.71  X  105 
N-m"3/2),  6.81  X  102  lb.-in._3/2  (7.49  X  105  N-nT3/2)  for  the  same 
respective  a/w  ratios.  Figures  4.4,  4.5,  and  4.6,  show  graphs  of  Kj„  vs. 
a/w. 

c)  Energy  Release  Rates 
Energy  release  rates  were  calculated  three  ways,  each  of  which  is 
described  below. 

1.  The  first  method  calculated  the  energy  release  rate  using 
the  corresponding  stress  intensity  value.  The  equations  used  are  as 
follows: 

Gj  =  1  '  v2  Kj2  For  unextended  a/w (4.1) 

l    2    2 
GIC  =      Kic  For  extended  a/w (4.2) 

These  values  are  listed  in  Tables  4.1,  4.2,  and  4.3. 

2.  The  second  method  calculated  the  energy  release  rate  by 
dividing  the  total  energy  to  failure  by  the  remaining  uncracked  area.  The 
total  energy  to  failure  was  found  by  calculating  the  area  under  the  load 
versus  LPD  trace  up  to  the  point  of  unstable  crack  growth.  This  point  was 
chosen  because  the  beam  was  failed  under  load  control.  The  uncracked  area 
was  calculated  as  follows: 

A  =  bw(l  -  a/w)(1.15) (4.3) 
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The  factor  1.15  was  used  to  compensate  for  the  rough  surface  (3).  The 
energy  release  rates  were  calculated  for  unextended  and  extended  a/w. 
Average  values  of  these  energy  release  rates  based  on  unextended  crack 
length  are  .248  lb. /in.  (43.45  N/m),  .196  lb. /in.  (34.34  N/m),  and  .155 
lb. /in.  (27.16  N/m)  for  approximate  initial  a/w  ratios  of  0.3,  0.5,  0.7. 
Values  based  on  extended  crack  length  by  CMOD  are  .308  lb. /in.  (53.96 
N/m),  .232  lb. /in.  (40.65  N/m),  .232  lb. /in.  (40.65  N/m)  for  the  same 
respective  a/w  ratios.  Values  based  on  extended  crack  lengths  by  IPD  are 
.280  lb. /in.  (49.06  N/m),  .207  lb. /in.  (36.27  N/m),  .183  lb. /in.  (32.06 
N/m)  for  the  same  respective  a/w  ratios.  Graphs  of  St*  calculated  this 
way  versus  a/w  are  shown  in  figures  4.7,  4.8,  and  4.9.  All  values  are 
listed  in  Table  4.4. 

3.  The  third  method  calculated  the  energy  release  rate  using 
the  J-integral  approach  in  which  the  slope  of  a  graph  of  the  energy  up  to 
the  point  of  instability,  U,  versus  a/w  is  used.  The  slope  of  this  line 
divided  by  the  effective  total  area,  that  is,  1.15  X  b  X  w,  gives  the 
energy  release  rate.  Plots  were  made  using  both  the  unextended  value  and 
the  extended  values  of  a/w.  These  are  shown  in  Figure  4.10,  4.11,  and 
4.12.  The  energy  release  rate  was  calculated  as  follows: 

Gi  ■  slope   For  unextended  a/w (4.4) 

d"w(1.15) 

GTr  =  s1°Pe   For  extended  a/w (4.5) 

IC   dw(1.15) 

Using  beam  with  a/w  of  0.5  and  0.7,  the  Gj  values  were  found  to  be  .251 
lb. /in.  (66.96  N/m)  for  unextended  a/w  values,  .207  lb. /in.  (55.22  N/m) 
for  extended  a/w  values  based  on  CMOD,  and  .222  lb. /in.  (59.22  N/m)  for 
extended  a/w  values  based  on  LPD.  Average  values  of  .229  (for  CMOD)  and 
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.251  (for  LPD)  are  used  in  comparing  to  the  other  methods.  It  is  seen  that 
these  values  are  fairly  consistant  for  nominal  a/w  values  of  0.5  and  0.7. 

A  summary  of  the  average  KjC  and  GjC  values  is  shown  in  Table  4.5. 
(d)  CTOD 

The  CTOD  was  calculated  using  an  equation,  based  on  William's  Stress 
Function,  which  was  developed  by  Go  and  Swartz  in  Reference  3,  which 
relates  the  CMOD  to  the  CTOD.  The  equation  is  as  follows: 

CTOD  =   R(w-a)   CMOD (2.3) 

R(w-a)  +  a  +  z 

R  =  Rotation  factor  taken  as  0.45  for  three  point  bending 
w  ■  beam  depth 
a  =  crack  depth 

z  =  .125  in.  (3.1  mm)  ■  knife  edge  thickness. 
CMOD  =  Crack  mouth  opening  displacement  at  unstable 
crack  growth 

The  CTOD  was  calculated  for  both  the  unextended  and  the  extended  crack 
lengths,  and  these  values  are  listed  in  Table  4.6. 

The  average  CTOD  values  based  on  unextended  crack  length  a/w  ratios 
of  0.3,  0.5,  0.7  are  8.92  X  10"4  in.  (2.27  X  10~2nm),  7.65  X  10"4  in  (1.94 
X  10"2mm),  6.75  X  10"4  in.  (1.72  X  10"2mm)  respectively.  For  extended 
crack  length  a/w  ratios  based  on  CMOD,  the  average  CTOD  values  are  6.41  X 
10"4  in.  (1.63  X  10"2mm),  6.18  X  10"4  in.  (1.57  X  10~2mm),  4.73  X  10-4 
in.  (1.20  X  10"2mm).  For  extended  crack  length  a/w  ratios  based  on  LPD, 
the  average  CTOD  values  are  7.49  X  10"4  in.  (1.90  X  10"2nm),  7.47  X  10"4 
(1.90  X  10"2mm),  5.63  X  10"4  in.  (1.43  X  10"2mm)  for  the  respective 
approximate  a/w  ratios  of  0.3,  0.5,  and  0.7.  The  CTOD  values  based  on 
extended  crack  length  give  consistant  results  regardless  of  the  initial 
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which  may  be  significant  in  determing  a  suitable  fracture  toughness 
parameter. 

4.3  Teflon  Beams 

The  six  teflon  beams  were  tested  as  previously  described.  The 
extended  a/w  was  found  and  the  stress  intensity,  energy  release  rate,  and 
CTOD  values  were  calculated  in  the  same  manner  as  in  the  precracked  beam. 
These  results  are  listed  in  Tables  4.7  and  4.8.  Average  values  are  listed 
in  Table  4.9. 

4.4  Strain  Beam 

The  beam  with  strain  gages  was  tested  as  previously  described.  The 
data  were  used  to  obtain  profiles  which  are  shown  in  Appendix  II.  From 
these  profiles  the  strain  reversal  point  was  located,  and  the  rotation 
factor,  R,  was  calculated  for  these  profiles  and  is  noted  on  each  profile 
in  Appendix  III.  The  values  found  proved  inconsistant,  so  the  value  of  R= 
0.45  as  found  by  Go  and  Swartz  in  Reference  3  was  used  in  the  calculation 
of  all  CTOD  values. 
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Table  4.1   Kj  and  (ij  Values  for  Unextended  a/w  in  Precracked  Beams 


Beam 
No. 

P-max 
(lb.) 

a/w 

unextended 

Kj  (lb.-- 

in, 

.  "3/2) 

fij  (lb. /in.) 

B-9* 

915 

.276 

7.90 

X 

10* 

.112 

B-10* 

905 

.330 

9.06 

X 

102 

.148 

B-ll 

955 

.307 

8.98 

X 

102 

.145 

C-l 

890 

.314 

8.53 

X 

102 

.131 

C-2** 

770 

.326 

7.62 

X 

102 

.105 

B-14 

480 

.506 

7.78 

X 

102 

.109 

B-16 

525 

.514 

8.71 

X 

10* 

.137 

B-17** 

390 

.521 

6.60 

X 

102 

.079 

C-3 

470 

.504 

7.58 

X 

102 

.103 

C-4 

445 

.525 

7.62 

X 

102 

.105 

B-18 

225 

.679 

6.38 

X 

10* 

.073 

B-19** 

190 

.685 

5.51 

X 

102 

.055 

B-20 

220 

.671 

6.05 

X 

10* 

.066 

C-5 

155 

.706 

4.89 

X 

10* 

.043 

C-6 

190 

.673 

5.27 

X 

10* 

.050 

Beams  loaded  at  a  faster  rate. 

** 
Beams  precracked  one  day,  and  failed  the  next. 

1  lb.  =  4.45  N 

1   lb.-in._3/2  =  3.48  X  10~2  N-mm~3/2 

1   lb. /in.   ■  0.175  N/mm 
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Table  4.2   KjC  an(j  GIC  Values  for  Extended  a/w  in 
Precracked  Beams  Based  on  CMOD 


Beam 

P-max 

a/w 

No. 

(lb.) 

extended 

B-9* 

915 

.450 

B-10* 

905 

.405 

B-ll 

955 

.430 

C-l 

390 

.480 

C-2** 

770 

.450 

B-14 

480 

.585 

B-16 

525 

.588 

B-17** 

390 

.618 

C-3 

470 

.560 

C-4 

445 

.575 

B-18 

225 

.773 

B-19* 

190 

.790 

B-20 

220 

.817 

C-5 

155 

.766 

C-6 

190 

.780 

KIC  (lb. -in.  "3/2)     aIC  (lb. /in.) 


12.68  X  102  .290 

11.09  X  102  .222 

12.53  X  102  .283 

13.25  X  102  .324 

10.67  X  102  .205 

9.85  X  102  .175 

10.88  X  102  .213 

8.68  X  102  .136 

8.93  X  102  .144 

8.85  X  102  .141 

9.65  X  102  .168 

8.93  X  102  .144 

10.23  X  102  .189 

6.62  X  102  .079 

8.45  X  102  .129 


Beams  loaded  at  a  faster  rate. 

** 
Beams  precracked  one  day,  and  failed  the  next. 

1  lb.  =  4.45  N 

1  lb.-in."3/2  =  3.48  X  10"2  N-mnf3/2 

1  lb. /in.  =  0.175  N/mm 


49 
Tab1e  4-3  Kic  and  GIC  Values  for  Extended  a/w  in 


Precracked 

Beams  Based 

on  LPD 

Beam 
No. 

P-Max 
(lb.) 

a/w 
Extended 

KIC  (lb. -in.  ~3/2) 
9.16  X  102 

Src  (lb. /in.) 

B-9* 

915 

.330 

.151 

B-10* 

905 

.360 

9.82  X 

lO2 

.174 

B-ll 

955 

.420 

12.19  X 

102 

.268 

C-l 

890 

.415 

11.21  X 

10* 

.226 

C-2** 

770 

.390 

9.06  X 

10* 

.148 

B-14 

480 

.540 

8.59  X 

10* 

.133 

8-16 

525 

.490 

8.13  X 

102 

.119 

B-17** 

390 

.622 

9.03  X 

102 

.148 

C-3 

470 

.420 

6.00  X 

102 

.065 

C-4 

445 

.535 

7.85  X 

102 

.111 

B-18 

225 

.730 

7.86  X 

102 

.111 

B-19* 

190 

.735 

6.78  X 

102 

.083 

B-20 

220 

.690 

6.51  X 

102 

.076 

C-5 

155 

.725 

5.30  X 

102 

.050 

C-6 

190 

.760 

7.63  X 

102 

.105 

* 
Beams 

loaded  at  a 

faster  rate. 

** 

Beams  precracked  one  day,  and  failed  the  next. 
1  lb.  =  4.45  N 

1  lb. -in. "3/2  =  3.48  x  10"2  N-nm"3/2 
1  lb. /in.  =  0.175  N/mm 
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Table  4.4   U  and  GjC  values  in  Precracked  Beams 


Beam 
No. 

U 
(lb. -in.) 

GIC  (lb/in.) 
(unextended  a/w) 

GIC  (lb/in.) 
(extended  CM00  a/w) 

SIC  (lb/in.) 
(extended  LPD  a/w) 

B-9* 

2.18 

.218 

.287 

.236 

B-10* 

1.96 

.212 

.239 

.222 

B-ll 

3.12 

.326 

.397 

.390 

C-l 

2.77 

.293 

.386 

.343 

C-2** 

1.76 

.189 

.232 

.209 

B-14 

1.51 

.222 

.264 

.238 

8-16 

1.44 

.215 

.253 

.205 

B-17** 

1.51 

.229 

.286 

.289 

C-3 

0.86 

.126 

.142 

.107 

C-4 

1.25 

.191 

.213 

.195 

B-18 

0.89 

.201 

.284 

.239 

B-19** 

0.61 

.140 

.210 

.167 

B-20 

0.73 

.160 

.289 

.171 

C-5 

0.46 

.114 

.142 

.121 

C-6 

0.72 

.160 

.237 

.217 

TT 

Beams  loaded  at  a  faster  rate. 
Beams  precracked  one  day,  and  failed  the  next. 
1  lb. -in.  =  113.03  N-nrn 
1  lb. /in.  =  0.175  N/mm 
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Table  4.5  Summary  of  Average  Values 


Nominal 
a/w 

Average 
Initial  a/w 

Average 

Final  a/w 

(CMOD) 

a/w 

% 
Chanqe 

Average 

Final  a/w 

(LPD) 

a/w 
% 

Change 

Group  1 

0.3 

0.311 

0.433 

42% 

0.383 

23% 

Group  2 

0.5 

0.514 

0.585 

14% 

0.521 

1% 

Group  3 

0.7 

0.683 

0.785 
CMOD 

15% 

0.728 

7% 

Average 


*IC 


Final  a/w   (lb. -in.  '3/2)  From  K 


GIC  (lb. /in.)   GIC  (lb. /in.)   GIC  (lb. /in. 


IC 


From  U 


From  J 


IC 


0.433 

1207 

0.265 

0.308 

0.585 

944 

0.162 

0.232 

.207 

0.785 

878 

0.142 
LPD 

0.232 

.207 

Average 

KIC 

GIC  (lb. /in.) 

GIC  (lb. /in.) 

GIC  (lb. /in.)* 

Final  a/w 

(lb. 

-in. 

iU) 

From  KjC 

From  U 

From  Jjq 

0.383 

1040 

0.193 

0.280 

-- 

0.521 

792 

0.115 

0.207 

0.251 

0.728 

681 

0.085 
Unextended 

a/w 

0.183 

0.251 

Average 


'IC 


Initial  a/w   (lb. -in. 


■3/2^ 


GIC  (lb. /in.)   GIC  (lb. /in.)   GIC  (lb. /in.) 


From  K 


IC 


From  U 


From  Jjq 


0.311 

841 

.128 

.248 



0.514 

766 

.106 

.196 

.251 

0.683 

562 

.057 

.155 

.251 

Based  on  initial  a/w  of  0.5  and  0.7. 
1  lb.-in."3/2  ■  3.48  X  10-2  N-mnf3/2 
1  lb. /in.  =  0.175  N/mm 
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Table  4.6   CMOD  and  CTOD  Values  for  Precracked  Beams 


Beam 
No. 

CMOD  (in.)  at 
instability 

1.99  X  IO-3 

CTOD  (in.)  for 
a/w  unextended 

10.24  X  IO-4 

CTOD  (in.)  for 

a/w  extended 

(CMOD) 

CTOD  (in.)  for 

a/w  extended 

(LPD) 

B-9* 

6.76  X  IO-4 

9.05  X  IO-4 

B-10* 

1.47  X  IO-3 

6.69  X 

lO"4 

5.59  X  IO"4 

6.23  X  IO"4 

B-ll 

1.80  X  IO*3 

8.64  X 

IO"4 

6.43  X  IO"4 

6.60  X  IO-4 

C-l 

2.30  X  IO"3 

10.86  X  IO"4 

7.22  X  IO"4 

8.53  X  IO-4 

C-2** 

1.78  X  IO"3 

8.17  X 

lO"4 

6.05  X  IO"4 

7.02  X  IO"4 

B-14 

2.70  X  IO"3 

7.90  X 

lO"4 

6.28  X  IO-4 

7.18  X  IO"4 

8-16 

2.92  X  IO"3 

8.36  X 

io-4 

6.73  X  IO-4 

8.93  X  IO"4 

B-17** 

3.02  X  10-3 

8.48  X 

io-4 

6.32  X  IO"4 

6.24  X  IO"4 

C-3 

2.33  X  10-3 

6.86  X 

IO"4 

5.85  X  IO"4 

8.54  X  IO"4 

C-4 

2.39  X  IO"3 

6.64  X 

IO"4 

5.73  X  IO"4 

6.45  X  IO'4 

B-18 

4.34  X  IO"3 

7.34  X 

IO"4 

5.52  X  IO"4 

5.97  X  IO"4 

B-19** 

3.80  X  lO'3 

6.28  X 

IO"4 

4.54  X  IO"4 

5.12  X  IO"4 

B-20 

4.38  X  IO"3 

7.63  X 

IO"4 

4.96  X  IO'4 

7.10  X  IO"4 

C-5 

3.60  X  10-3 

5.48  X 

IO"4 

4.20  X  IO*4 

5.06  X  IO'4 

C-6 

4.08  X  IO"3 

7.05  X 

IO"4 

4.44  X  IO-4 

4.90  X  IO"4 

AVERAGE 

7.77  X 

IO"4 

5.77  X  IO"4 

6.86  X  IO"4 

** 


Beams  loaded  at  a  faster  rate. 
Beams  precracked  one  day,  and  failed  the  next. 
1  in.  =  25.40  mm 
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Table  4.7      Kj  an(j  Gj  Values  in  Teflon  Beams 


B< 
N( 

iam 
). 

•15 

P-Max 
(lb.) 

570 

a/w 

Unextended 

Kj   (lb.-- 
5.16 

in.  -3'2) 
X  102 

«!   (lb. /in.) 

C- 

.293 

.048 

C- 

■16 

570 

.325 

5.63 

X  102 

.057 

C- 

■17 

385 

.455 

5.41 

X  102 

.053 

C- 

■18 

290 

.514 

4.81 

X  102 

.042 

c- 

■19 

160 

.625 

3.74 

X  102 

.025 

c- 

Be 

•20 

i. 

•15 

95 

P-Max 
(lb.) 

570 

a/w 

.670 

1  Extended 
CMOD 

2.60 
KJC  (lb.- 

9.35 

X  102 
■in.  -3/2) 

X  102 

.012 
GIC  (lb. /in.) 

C- 

.510 

.157 

C- 

•16 

570 

.513 

9.43 

X  102 

.160 

c- 

•17 

385 

.617 

8.76 

X  102 

.138 

c- 

■18 

290 

.670 

7.95 

X  102 

.114 

c- 

■19 

160 

.825 

9.32 

X  102 

.157 

c- 

Bf 

N( 

■20 
:am 
). 

■15 

95 
P-Max 
(lb.) 

570 

a/vt 

.910 
Extended 

(LPD) 

12.39 
KIC   (lb.- 

13.32 

X  102 
•in.   -3/2) 

!  X  102 

.276 
GIC  (lb. /in.) 

c- 

.625 

.320 

c- 

■16 

570 

.693 

17. 01 

i  X  102 

.524 

c- 

•17 

385 

.717 

12.71 

.  X  102 

.291 

c- 

•18 

290 

.732 

10.22 

:  x  io2 

.188 

c- 

■19 

160 

.805 

8.21 

.  X  io2 

.121 

c- 

•20 

95 

__ 

,. 

__ 

1   lb.   =  4.45  N 

1   lb.-in._3/2  -  3.48  X  IO"2  N-rmf3/2 

1   lb. /in.   =  0.175  N/mm 
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Table  4.8   U,  GIC>  CMOD  and  CTOD  Values  for  Teflon  Beams 


Beam 
No. 

U       GIC  (lb/in)      GIC  (lb/in) 
(lb. -in.)   Unextended  a/w   Extended  CMOD  a/w 

GIC  (lb/in) 
(Extended  LPD  a/w) 

C-15 

2.40 

.246 

.355 

.464 

C-16 

3.30 

.354 

.491 

.778 

C-17 

1.75 

.232 

.331 

.448 

C-18 

1.18 

.176 

.259 

.319 

C-19 

.86 

.167 

.358 

.321 

C-20 

.39 

.084 

.309 

— 

Beam 
No. 

CMOD  (in.)  at 
instability 

1.86  X  IO"3 

CTOO  (in.)  for 
a/w  unextended 

9.21.X  IO"4 

CTOD  (in.)  for 

a/w  extended 

(CMOD) 

CTOD  (in.)  for 

a/w  extended 

(LPD) 

C-15 

5.38 

X  IO"4 

3.81  X  IO"4 

C-16 

1.92  X  KT3 

8.84  X 

10-4 

5.51 

X  IO"4 

3.08  X  IO"4 

C-17 

2.84  X  IO"3 

9.52  X 

io-4 

5.96 

X  IO"4 

4.13  X  IO"4 

C-18 

3.24  X  IO"3 

9.28  X 

io-4 

5.66 

X  IO"4 

4.42  X  IO"4 

C-19 

4.40  X  KT3 

9.00  X 

io-4 

3.71 

X  IO"4 

4.18  X  IO"4 

C-20 

4.38  X  KT3 

7.65  X 

IO"4 

1.81 

X  IO"4 

AVERAGE 

5.11  X 

IO"4 

4.67 

X  IO"4 

3.92  X  IO"4 

1  lb. 

-in.  =  113.03  N-m 

1  lb. 

/in.  =  0.175  N/mm 

1  in. 

=  25.40mm 
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Table  4.9  Summary  of  Average  Values  of  KTr  and  Grr  for  Teflon  Beams 


IC 


JIC 


Nominal 
a/w 

Average 
Initial  a/w 

0.309 

Average 
Final  a/w 

(CMOD) 

a/w 

% 
Chanqe 

66% 

Average 

Final  a/w 

(LPD) 

a/w 

% 

Chanqe 

0.3 

0.512 

0.659 

113% 

0.5 

0.485 

0.644 

33% 

0.725 

49% 

0.7 

0.648 

0.868 
CMOD 

59% 

0.805 

24% 

Average 
Final  a/w 

(lb.-in.-3/2) 

S] 

:c  (lb. /in.) 
From  KyC 

GIC 

(lb. /in.) 
From  U 

0.512 

9.39  X  102 

0.159 

0.423 

0.644 

8.36  X  102 

0.126 

0.295 

0.868 

10.86  X  102 

LPD 

0.217 

0.334 

Average 
Final  a/w 

KIC 
(lb.-in.-3/2) 

SIC  (lb. /in.) 
From  K.p 

GIC 

(lb. /in.) 
From  U 

0.659 

15.19  X  102 

0.422 

0.621 

0.725 

11.47  X  102 

0.240 

0.384 

0.805 

8.21  X  102 

Unextended 

0.121 
.  a/w 

0.321 

Average 
Initial  a/w 

KIC 
(lb. -in."3'2) 

GIC  (lb./in."3/2) 
From  KjC 

GIC 

(lb. /in.) 
From  U 

0.309 

5.39  X  102 

0.053 

0.300 

0.485 

5.11  X  102 

0.047 

0.204 

0.648 

3.17  X  102 

0.019 

0.126 

1   lb.-in.-3/2  =  3.48  X  10"2  N-imf 
1   lb. /in.   ■  0.175  N/mm 


3/2 


56 
Chapter  5 

SUMMARY  AND  CONCLUSIONS 

The  experimental  results  are  interpreted  and  summarized  as  follows: 

1.  The  extended  crack  lengths  in  precracked  beams  based  on  the  CMOD 
compliance  calibration  curve  were  found  to  be  39%,  145t  and  15%  higher 
for  nominal  a/w  ratios  of  0.3,  0.5,  and  0.7  respectively  than  the 
unextended  crack  lengths.  For  those  based  on  the  LPD  curve,  the 
extended  crack  lengths  were  found  to  be  23%,  1%  and  7%  higher  for  the 
same  nominal  a/w  ratios.  The  teflon  beams  exhibited  larger  variations 
of  56%,  33%,  and  34%  based  on  the  CMOD  curve  and  113%,  49%,  and  24% 
based  on  the  LPD  curve  for  nominal  a/w  ratios  of  0.3,  0.5,  and  0.7 
respectively. 

2.  The  Kj  vaiues  for  precracked  beams  were  consistantly  higher  than  for 
teflon  beams,  supporting  earlier  reports.  The  Kj  values  for  extended 
crack  lengths  based  on  CMOD  were  44%,  23%,  and  56%  higher  than  the 
unextended  crack  lengths  for  nominal  a/w  ratios  of  0.3,  0.5,  and  0.7. 
Based  on  LPD,  the  crack  lengths  were  22%, 3%,  and  21%  higher.  Although 

tne  Kj  value  is  higher  for  smaller  a/w  ratios,  it  is  fairly  consistant 
for  a/w  ratios  of  approximately  0.5  and  0.7. 

3.  The  energy  release  rates  showed  the  following  trends: 

(a)  Those  found  using  stress  intensity  values  were  much  lower  than 
those  found  by  the  energy  method  for  both  the  unextended  crack 
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length  and  the  extended  crack  length  based  on  CMOD.   The 
values  found  based  on  LPD  were  only  slightly  lower. 

(b)  Based  on  the  data  from  nominal  a/w  ratios  of  0.5  and  0.7  there 
is  good  agreement  between  the  energy  release  rates  based  on 
the  energy  approach  and  the  J-integral  approach  for  unextended 
a/w  values  and  extended  a/w  values  based  on  CMOD  as  well  as 
LPD. 

(c)  Energy  release  rates  found  using  extended  a/w  ratios  based  on 
CMOD  tended  to  be  higher  than  those  based  on  LPD  except  those 
found  by  the  J-integral  approach  which  only  differed 
slightly. 

4.  The  CTOD  values  found  using  the  extended  crack  lengths  based  on  CMOD 
and  LPD  gave  very  consistant  results  °f  5.77  X  10"4  in.  (1.47  X  10"2mm) 
and  6.86  X  10"4  in.  (1.74  X  10"2mm)  respectively  with  less  than  10* 
average  deviation.  The  CTOD  based  on  the  unextended  a/w  values  gave  an 
average  value  of  7.77  X  10"4  with  less  than  15%  average  deviation.  The 
CTOD  values  found  for  teflon  beams  were  only  slightly  less.  The 
consistancy  of  the  CTOD  demonstrates  that  it  may  be  an  excellent 
indicator  of  fracture  toughness. 

Recommendations  for  future  work  include  the  investigation  of  size  variation 
of  the  speciman  and  continued  statistical  evaluation  of  the  methods 
presented  here. 
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APPENDIX  II 

Sand  and  Aggregate  Properties 


A  specific  gravity  test  was  performed  on  the  sand  as  per  ASTM  C128-79. 
Briefly,  the  following  procedure  was  used: 

Weigh  approximately  500  grams  of  oven  dry  sand  and  place  into  a 
pycnoneter.  Add  distilled  water  to  fill  the  pycnometer  to  about 
75%  to  90%  capacity.  Using  a  vacuum  and  a  burner,  boil  the  air 
out  of  the  sample  and  water,  shaking  the  bottle  continuously. 
(This  takes  about  20  minutes.)  Cool  to  room  temperature  and  add 
de-aired,  distilled  water  to  the  calibrated  mark  on  the 
pycnometer.  Weigh  the  bottle,  water  and  sand,  and  measure  the 
temperature  of  the  sample.  Empty  the  contents  of  the  pycnometer 
into  an  evaporating  dish  using  distilled  water  to  remove  all  of 
the  grains  from  the  bottle.  Fill  the  pycnometer  with  distilled 
water  to  the  calibrated  mark  and  weigh  (Make  sure  the  temperature 
is  the  same  as  recorded  before).  Place  the  sand  and  water  in  a 
dryinq  oven  until  all  the  moisture  has  been  removed,  about  24 
hours.  Weigh  the  dried  sand.  Calculate  the  specific  gravity  using 
the  following  equation: 

S.G.  (at  temp,  measured)* wt.  dry  sand 

(wt.flask+water)+(wt.dry  sand )-(wt. flask, sand, water) 

The  test  was  performed  on  two  samples.  Sample  1  had  a  S.G.  of  2.658  and 
smple  2  had  a  S.G.  of  2.651.  The  average  of  these  two  was  used 
(S.G. =2. 655) 


A  specific  gravity  test  was  performed  on  the  aggregate  as  per  ASTM 
C127-81.  Briefly,  the  following  procedure  was  used: 

For  3/4"  aggregate,  soak  6.6  pounds  of  aggregate  in  water  for  24 
hours  prior  to  testing.  After  soaking,  remove  the  sample  from 
the  water  and  dry  the  surface  of  the  particles.  Weigh  the 
saturated  surface-dry  sample.  Put  the  sample  in  a  wire  bucket 
with  a  small  enough  mesh  that  none  of  the  aggregate  will  fall 
through.  Weigh  the  sample  submerged.  Remove  the  sample  and 
place  in  an  evaporating  dish  and  oven  dry  for  24  hours.   When 
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sample  is  completely  dry,  weigh  again.  Use  the  following  formula 
to  calculate  the  specific  gravity: 

S.G.»  (wt. oven  dry) 

(wt. sat., surf ace  dry)-(wt. sat., submerged) 

This  test  was  performed  on  two  samples.  Sample  1  had  a  S.G.  of  2.54  and 
sample  2  had  a  S.G.  of  2.57.  The  average  of  these  two  was 
used.  (S.G. =2. 56)  A  sieve  analysis  was  done  on  both  the  sand  and  the 
aggregate  as  per  ASTM  C136-82. 

For  the  sand,  500  gram  samples  were  used  for  each  test.  The 
ovendry  sample  was  placed  into  a  standard  sieve  set  consisting  of 
the  following  sieve  sizes:  .150mm,  .300mm,  .600mm,  1.18mm, 
2.38mm.  The  sieves  were  shaken  for  10  minutes  in  the  mechanical 
shaker,  and  the  amount  retained  on  each  sieve  was  recorded.  The 
fineness  modulus  was  calculated  to  be  2.91,  and  the  gradation 
curve  was  plotted. 

The  aggregate  sample  was  11  pounds  as  recommended  and  was  broken 
into  five  parts.  Each  of  the  five  partial  samples  was  put 
through  the  sieve  test  using  the  following  sieve  sizes: 
#10(2. 00mm),  #8(2. 38mm),  #4(4.76mm),  3/8"(9.52mm),  l/2"(12.70mm), 
3/4  (19.10mm).  As  with  the  sand,  they  were  shaken  with  the 
mechanical  shaker  for  10  minutes  and  the  gradation  curve  was 
plotted. 


In  both  the  sand  and  aggregate  sieve  analysis,  an  analysis  of  any  material 
finer  than  the  75X10"6m  sieve,  that  is  the  dust,  was  neglected. 

A  dry  rodded  unit  weight  was  performed  on  the  aggregate  as  per  ASTM 
C29-78. 

Find  the  volume  of  a  bucket  by  weighing  the  bucket  empty,  then 
weighing  it  full  of  water  at  room  temperature,  and  then  dividing 
the  weight  of  the  water  by  the  density  of  water  at  room 
temperature.  Fill  the  bucket  1/3  full  and  rod  25  times.  Add 
another  1/3  and  rod  25  times  penetrating  the  first  layer  by  1-2 
inches.  Repeat  with  the  third  layer.  Level  off  the  bucket  and 
weigh  the  bucket  and  aggregate.   Divide  the  weight  of  the 
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aggregate  by  the  volume  of  the  bucket  to  get  the  dry  rodded  unit 
weight. 

This  test  was  performed  two  times  and  the  average  dry  rodded  unit  weiqht 
was  94  lbs. /ft.3 
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Appendix  III 
STRAIN  PROFILES 
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Appendix  IV 
RAW  DATA 

Scale  for  all  raw  data  graphs  is  as  follows: 

X-axis    1  square  (as  noted  on  p. 75)  =  2  cm 
Y-axis     1  square  ■  2  cm 

Scale  factors  to  determine  actual  values  of  load  (Y-axis)  and  displacement 
(X-axis)  are  found  on  pages  13  and  14.  Each  graph  is  labeled  CMOD  or  LPD 
to  identify  what  kind  of  displacement  is  represented. 
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Cylinder  Test  Data 


Compressive  Strength  in  PS  I 


Batch  A: 


4074 

4711* 

7569 

7512 

7172 

8715 

8573 

8191 


Batch  B: 


8078 
7993 


9181 
7116 


7908 
8446 
8361 


8474 


Tensile  Strength  in  PSI 


Batch  A: 


601 


Batch  B: 


665 


Testing  machine  not  working  properly  so  these  values  are  disregarded. 


1.0  PSI  =  6.89  X  10   MPa 
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Abstract 

An  extensive  experimental  program  has  been  started  in  an  attempt  to 
verify  the  relationships  between  fracture  parameters  for  plain  concrete  in 
bending.  Results  of  a  pilot  program  indicated  promising  results  if  the 
crack  lengths  at  the  onset  of  unstable  crack  growth  were  considered  in  the 
evaluation  of  the  parameters.  It  was  also  demonstrated  that  the  crack  tip 
opening  displacement  (CTOO)  may  be  a  valid  fracture  criterion.  This  study 
attempts  to  verify  these  relationships  with  statistical  confidence  for 
three  sizes  of  beams.  The  program  presented  here  includes  only  one  size 
of  beam,  3  in.  wide,  4  in.  deep  with  a  15  in.  span.  Future  studies  will 
be  made  on  larger  beams. 

These  small  beams  were  precracked  using  a  modified  compliance 
calibration  technique  and  loaded  to  failure.  Plots  of  load  versus  crack 
mouth  opening  displacement  (CMOD)  and  load  versus  load  point  displacement 
(LPD)  were  obtained  simultaneously.  Five  beams  of  each  crack  depth  to 
beam  depth  ratio  (a/w)  of  0.3,  0.5,  and  0.7  were  tested  in  three  point 
bending.  Two  teflon-insert  (prenotched)  beams  were  loaded  to  failure. 
Additionally,  one  beam  was  instrumented  with  strain  gages  in  order  to 
obtain  strain  profiles  ahead  of  the  crack  tip. 

Results  presented  include  stress  intensity  values  (KI))  energy  release 
rates  (Gj)  t>aseC|  on  three  methods  of  analysis,  and  CTOD  values. 
Comparisons  are  made  between  these  values  based  on  unextended  crack 
lengths  and  on  extended  crack  lengths. 
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